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A SPACE -TIME TENSOR FORMULATION FOR 
CONTINUUM MECHANICS IN GENERAL CURVILINEAR, MOVING, 
AND DEFORMING COORDINATE SYSTEMS 


Lee M. Avis 

Langley Research Center 
SUMMARY 

Tensor methods are used to express the continuum equations of motion in general 
curvilinear, moving, and deforming coordinate systems. The space-time tensor formu- 
lation is applicable to situations in which, for example, the boundaries move and deform. 
Placing a coordinate surface on such a boundary simplifies the boundary -condition treat- 
ment. The space -time tensor formulation is also applicable to coordinate systems with 
coordinate surfaces defined as surfaces of constant pressure, density, temperature, or 
any other scalar continuum field function. The vanishing of the function gradient compo- 
nents along the coordinate surfaces may simplify the set of governing equations. In 
numerical integration of the equations of motion, the freedom of motion of the coordinate 
surfaces provides a potential for enhanced resolution of the continuum field function. The 
space -time tensor formulation has been applied to numerical simulations of the atmo- 
sphere and the oceans. 

Two tensor expressions of inertial velocity which provide a convenient flexibility in 
formulating equations of motion are presented. One inertial velocity tensor, in contra- 
variant form, has space components equal to the (nontensor) velocity relative to the coor- 
dinate system and the time component equal to an arbitrary constant (the metric tensor 
providing the velocity of the coordinate system). The other inertial velocity tensor, in 
contravariant form, has a vanishing time component. Also, a tensor of coordinate system 
velocity relative to inertial space exists in the space -time. Expressions for the covariant, 
absolute, and comoving derivatives of general space-time tensors are derived in terms of 
the spatial Christoff el symbols and the coordinate system inertial velocity. Expressions 
for the rates of change of the elements of the space -time metric tensor are presented. 

The equations of motion of a material continuum in the space -time tensor formalism are 
derived from generalized conservation principles from both Lagrangian and Eulerian 
viewpoints. An example problem of an incompressible, inviscid fluid with a top free sur- 
face is considered, where the surfaces of constant pressure (including the top free surface) 
are coordinate surfaces. 



INTRODUCTION 


This paper presents a space -time tensor formulation for nonrelativistic continuum 
mechanics which is applicable to general time -dependent geometry and coordinate systems 
moving with respect to each other. A space-time metric tensor which preserves the 
invariance of the nonrelativistic spatial distance between neighboring simultaneous events 
is defined. Time assumes the absolute character appropriate to Newtonian mechanics. 

The space -time tensor formulation has been applied (ref. 1) to numerical simulation 
of atmosphere and ocean dynamics employing a coordinate system rotating with the Earth 
and having quasi-horizontal coordinate surfaces variable in space and time. This formu- 
lation, by virtue of its generality and compactness, has been found to be a useful device. 

In general, problems which might profitably be expressed in the space-time tensor 
formalism include those with moving and deforming boundaries and/or those in which one 
may choose as coordinate surfaces, surfaces of constant pressure, temperature, or other 
physical parameters. In the former case, the treatment of boundary conditions is facili- 
tated; and the latter may lead to material simplification of the governing equations. Also, 
tracking the chosen physical parameter with coordinate surfaces is potentially a means 
of improving the resolution of the parameter in numerical integration of the equations of 
motion, without decreasing the space and time numerical step sizes. Some examples of 
suitable coordinate -surface -defining functions are density, temperature, and pressure in 
simulations of the troposphere; and pressure (in many cases) in simulations of a general 
fluid in a gravitational field. 
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SYMBOLS 

determinant of matrix of Aj^j 

covariant and contravariant elements of space -time metric tensors 

covariant and contravariant elements of space -time metric tensors for a 
Riemann space -time which is a subspace of Euclidean space-time 

_3 

z component of fluid comoving acceleration 
= ^N+1, N+1 " ^ 

source term tensor, defined by equation (66) 
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DO 


absolute differential 

space -time tensor of force per imit volume 




Q!/3 




a/3 

,g 




K 

N 

N 

^i 

O 


covariant and contravariant elements of N-space metric tensor 


■,g P covariant and contravariant elements of a subspace of Euclidean 3-space 


determinants of 


and 


determinant of matrix of covariant elements of metric tensor in subspace 
of coordinates 

velocity of x coordinate system relative to inertial z system, 

■(*)., • ‘1 > 

constant between zero and unity 

number of spatial dimensions in Riemann space-time, N = 1,2,3 
number of dimensions in subspace of Euclidean 3-space, N = 1,2 


covariant vector normal to bottom boundary 


origin of Riemann coordinates 


pressure 




pressure at top free surface, atmospheric pressure 


inertial velocity of coordinate system, S/ > = - 


(y) 


9t 1 0 


0 component of inertial velocity of coordinate system in subspace of 
Euclidean 3 -space 



\ ds aloi^ geodesic from point O to point 
= yAjj dy* dy' 

= dZa 

general space -time tensor 
time 

inertial velocity, 

inertial velocity, dy^dt 

(u^ - J^)n- at bottom boundary 

inertial velocity of a particular coordinate system, the y system of 
appendix B; defined in appendix D (eqs. (D 9 ) and (DIO)) 

Riemann coordinates 

spherical polar coordinates fixed with respect to the rotating Earth; 
12 S 

x^, X , and x are south, east, and radially outward, respectively, 
and dx^^ = X dt 

general curvilinear, moving and deforming coordinates 
any coordinates such that -A-^j(2;) ” ®j(z) locally 

inertial Cartesian coordinate system over Euclidean space -time of 
four dimensions 

geocentric Cartesian coordinate system (approximately inertial) 


Kronecker delta, equal to zero for i 7^ j and equal to unity for i = j 



e 

X 

p 

S2 

Ct) 


U) 


ij 


gravitational potential of Earth 
arbitrary nonvanishing constant 
variable of integration 

mass density 

angular velocity of Earth’s rotation relative to inertial space 

angular velocity of rotation of pressure gradient relative to z coordinate 
system 


space -time vorticity, 



Notation: 



Christoffel symbol of the second kind in space-time 


Christoff el symbol of the second kind in N-space 


Christoffel symbol of the second kind in Riemann space -time which is 
a subspace of Euclidean space -time 


( ) 


() 


upper parenthetical term is expressed in coordinate system labeled by 
lower parenthetical symbol (upper parentheses may be absent) 



y^ coordinates are held fixed for the partial differentiation 


bracketed expression is evaluated in N-space 
parenthetical expression is evaluated at point O 


covariant differentiation 



determinant of matrix 





defined as zero 


defined as zero 


Indices; 


The range of Roman indices 
bars, 1 to N. 

i 

to 

u 

is 

1 

to 

N + 1; of Greek indices without 

The range of Roman indices 
1 to N. 

i 

to 

u 

is 

1 

to 

N + 1; of Greek indices a to 

The range of Greek indices 

a 

to 

w 

is 

1 

to 

(N - N). 


The index b refers to the bottom boundary. 

The summation convention observed herein is that a lowercase letter occurring 
exactly twice in a term, once as a lower index and once as an upper index, denotes sum- 
mation over those indices. 


ANALYSIS 


Introductory Remarks 


The tensor formulation developed in the analysis applies to coordinate transforma- 
tions of the form 


y' = y'(y*, • • •. y“ y^^i 
y' = y^, 




( 1 ) 




which are transformations of the coordinates of an event in space-time, where N is the 
number of spatial dimensions and 


-N.l ^ yN.l ^ 


( 2 ) 
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The symbol X is an arbitrary nonvanishing constant and t is the time elapsed since 
an arbitrary reference time. The functions y^ and y^ are required to be single valued 
and continuous with continuous derivatives through the third order. Such transformations 
and coordinate systems are called "allowed" herein. 

A brief outline and description of the development sequence is given as follows: 

"The Space-Time Metric Tensor" 

A space-time interval in a local Cartesian coordinate system, which establishes a 
space-time metric tensor over a Riemann "space," is introduced. The full four- 
dimensional space-time is Euclidean, but the generality of Riemann space-time is retained 
for application to surfaces with Gaussian curvature. Appendix A develops the relationship 
of the Christoffel symbols in Riemann space-time to those in Euclidean space-time. 

"Inertial Velocity Tensors" 

Two tensor expressions of inertial velocity which provide a convenient flexibility in 
formulating equations of motion are presented. One inertial velocity tensor has the time 
component X and the N contravariant space components equal to the (nontensor) veloc- 
ity relative to the coordinate system. The other inertial velocity tensor, in contravariant 
form, has a vanishing time component. The difference of the latter and former tensors is 
the coordinate system velocity tensor (which, upon coordinate transformation, applies to 
the coordinate system transformed to). 

"Space-Time Christoffel Symbols" 

The Christoffel symbols for the Riemann space -time are expressed in terms of the 
coordinate system velocity and the Christoffel symbols for N-space. 

"The Space-Time Tensorial Derivatives" 

The covariant, absolute, and comoving derivatives of a general tensor are derived by 
use of the expressions for the space-time Christoffel symbols. The covariant and 
comoving derivatives are then specialized for first- and second-rank tensors. 

"The Rates of Change of the Metric Tensors" 

The partial time derivatives of the covariant and contravariant space -time metric 
tensors and N-space metric tensors are derived. 

"Equations of Motion of a Material Continuum" 

The equations of transport for a general continuum tensor field are derived from 
generalized conservation principles, and then specialized to various forms of the equations 
of motion. A space -time vorticity tensor is introduced, and the equations of motion are 
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expressed in terms of vorticity and kinetic energy. The Eulerian and Lagrangian time 
derivatives of vorticity are expressed in terms of the fields of velocity, vorticity, and 
acceleration. 

These results are applied in appendix B to an example problem, employing constant- 
pressure coordinate surfaces, of an ideal fluid moving over the Earth. Appendixes C 
and D are supplements to appendix B. 

The Space -Time Metric Tensor 

This section introduces a space -time metric tensor which renders invariant the 
spatial part of the space-time interval as well as the total space-time interval. 

2 

The square of the space -time interval (ds) is defined by 

(ds)^ = dz. dz^ (3) 

where the N coordinates z“ are of a local Cartesian system regarded as fixed in 

inertial space and z^'*'^ = Xt. Transforming the z^ coordinate increments to any 
allowed coordinate system y yields 


9z i 


(4) 


where the spatial components of the space -time metric tensor can be written 



dz 

y 


9y^ 9y^ 9y^ 


!!y ^ 

9y“ 9y^ 


(5) 


since 


dz 


N+1 




and, in general, 


9y“ 9y“ 


'( 6 ) 
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The other components of A. 


are 


mn 


■'^a.N+l ■ ^N+l,a 


az^ 

9y® 8y® \ 


(7) 


^N+1,N+1 


gyN+1 gyN+1 


1 = 





(8) 


where, it might be noted, 



is the Cartesian N-vector of the velocity of the 


y coordinate system relative to the inertial z system. Equation (5) states that the 
spatial components of the covariant space -time metric tensor are the components of the 
covariant N-space metric tensor, or 

The time coordinate increment transforms by invariance, from equation (2); thus, 
rearranging equation (3) yields 


dz^ dz>^ = (ds)2 - dz^^^ dz^+^ = (ds)2 - (\)2(dt)2 = (ds)2 


which is invariant. 


The y coordinate system velocity 



occurring in equations (7) and (8) is 


expressed in y coordinates by 



\ _ dz^ _ 

dz'y 

dy" 

_ dz^ _ 

9z^ 

9y^^ 

dz^ 

3z^ 

[9y“] 


1/3 


dt 

dt 

9y“ 

dz^ 

dt 

9y“ 



or 



(9) 
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where 


is the velocity of the y coordinate system relative to the z system, 


expressed in y coordinates. Let 




, is! 

’“<5 8y“ 8y^ 


Then, from equations (5) and (7) to (11), 


^a/3 ^Q!/3 


A — A — ^ o- 

■^0!,N+1 “ ■^N+1,0' “ X ^acj) ° 


^N+1,N+1 “ (x) ^ctl3 + 1 


The components of the contravariant space -time metric tensor are given by 


9Zk 9z^ 




^Q!,N+1 ^ ^N+1,0! ^ 


-1 s“ 


^N+1,N+1 ^ ^ 


where g ^ is the contravariant N-space metric tensor. 


10 



It should be noted that the coordinate system velocity given by equation (10) is 

not a space -time tensor, because the time component has not been defined. A space -time 
tensor of coordinate system velocity is introduced in the following section. 


Inertial Velocity Tensors 

The space -time tensor 


dt 


(14) 


where dy^ 
space -time 


is the coordinate increment tensor associated with two neighboring events in 

and dt = i dy^, transforms to the local Cartesian inertial z system as 
\ 


V 


j _ d^ 

(z) dt 


the spatial components of which constitute the Cartesian N-vector of inertial velocity. 

The time component X can be interpreted as a constant uniform velocity normal to all 
spatial directions. If the two events in equation (14) are the presence of a certain particle, 
then is an inertial space -time velocity of the particle. 

Another inertial velocity space -time tensor is defined by the contravariant vector 
transformation from the z coordinate system to any allowed system y of the Cartesian 
space-time vector 




U 


N+1 

(z) 


= 0 


Thus, 




dz^ 9y^ 
dt 92 ^ 



(15) 
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which has a vanishing time component. From equation (10), the spatial components 
of are 


u“ = V® + s“ 


which expresses inertial velocity as the sum of coordinate system velocity S® and 
velocity relative to the coordinate system. 


In equation (15), the term 



is the difference between two space -time 


tensors and is therefore a space -time tensor. The tensor of coordinate system inertial 
velocity is defined as 



(16) 


and equation (15) becomes 

(17) 

The relationship among the velocities U^, V^, and under coordinate transformation 
is illustrated in sketch (a), which represents a two-dimensional case. The curve FG is 



Sketch (a) 
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the path taken by a point of fixed spatial y coordinate as seen in the z coordinate sys- 
tem. The line Xc represents of a particle at point A. The line AE, drawn 

tangent to the curve at point A, intersects BC to determine the spatial component 

of the y coordinate system velocity, represented by BE. The point E partitions the 
spatial component of inertial velocity y coordinate system velocity bIi 

and the particle velocity relative to the y system, expressed by 





+ 



or, in y^ lanits, from eqviations (9) and (10), 


"(y) = ^(y) " ®(y) 


The time components of U, V, and S are similarly related, since ^ 


'^(y) = -"(y) = 


TT^ — 4- 

%) " '^(y) * (y) 


Hence, 




1 2 

holds for any allowed transformation from the coordinate system z ,z of sketch (a). 

Lowering the indices of and and applying equations (12) and (17) yields, 
for any allowed coordinate system, 




a 


a 



Un ^.1 = 

X 0 

1 





> 

(18) 

= 0 

®N+1 ^ 

-X 
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The existence of two inertial velocity tensors V and U provides a convenient 
flexibility in formulating the dynamical equations in curvilinear, moving coordinate sys- 
tems. The space components of constitute the velocity relative to the coordinate 
grid, which is a convenient choice of velocity representation to use in the convective 
terms of the dynamical equations of a material continuum, for example. On the other 
hand, an advantage of representing inertial velocity by U is that the time component 
of vanishes. Inertial velocities might be represented by U at some places and 
by V at other places, whichever is more eonvenient, in the problem formulation. A 
given velocity representation can be tailored to one’s needs by use of equation (17). 


Space-Time Christoffel Symbols 

In preparation for the discussion of covariant and comoving derivatives, the space- 


time Christoffel symbols of the second kind 




A 


of the N- space Christoffel symbols 
velocity. 

The Christoffel symbols of the second kind, defined as 


are derived in this section in terms 
and the coordinate system inertial 


f. ^ '] 

U ^ \ 3y^ 9y^ 9y^/ 


(19) 


transform from an allowed coordinate system y to any allowed system y by (p, 48, 
ref. 2) 



] 

p q 


9y 


9yP dy^ I 9y' 


(y) 9y^ 9y^ 9y*^ 3yJ 


9^yj 

9y™ 9y^ 


(20) 


The second expression on the right-hand side implies that 


1 

j k 


is not a tensor. In 


Euclidean "space," such as four -dimensional space-time, it is always possible to define 

N+1 

a coordinate system z such that (ds)^ = ^ (dz^ dz^) everywhere (p. 583, ref. 3). 


i=l 


The metric tensor components are then constant and the Christoffel symbols vanish in 
the z coordinate system. Therefore, in Euclidean space -time, from equation (20), 
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(21) 




_ 9y"' 


2d 




az^ aym 9yH 


Equation (21) can be generalized for application to non-Euclidean Riemann space -time by 
use of Riemann coordinates x^, which are defined by (p, 584, ref. 3) 

= sp^ (22) 

where p^ is the unit tangent vector dy^/ds at the origin (xi = o) of a geodesic joining 
the origin and the point yJ and s is the geodesic ’’distance" between the two points, or 

px^(y^) 

s = \ ds (along geodesic) 

V=0 


From reference 3 (p. 584), Riemann coordinates can be defined in any Riemann space at 
any given origin, and at the origin all Christoffel symbols vanish. Therefore, from 
equation (20), at the origin of the x system. 



(y) dxi 3y^ ay’’ 


(23) 


The vanishing of the Christoffel symbols in Riemann coordinates at the origin often is 
used to simplify mathematical arguments. 

It is shown in the following paragraph that dx^'*'^ = X dt, as required by equation (2) 
for the X coordinate system to be allowed, 

A geodesic is uniquely defined by any given point on the geodesic, the direction (unit 
tangent vector) dy’/ds of the geodesic at the given point, and the equation (p. 583, 
ref. 3) 


(ds)2 



^ 0 

ds ds 


or 


ds \ds / 


= 0 


(24) 


(25) 
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where D is the absolute differential operator. In Riemann coordinates, equation (24) 
becomes, by equation (22), 


dV 

(ds)2 


= 0 


(26) 


The non-Euclidean space -time of interest here is a subspace of Euclidean space-time. 
It is shown in equations (A19) to (A21) of appendix A that 



= 0 


(27) 



(28) 


in any allowed coordinate system over Riemann space-time. From equations (24) and (27), 


.2 N +1 . 2 , 

=0 


(ds)' 


(ds)' 


(29) 


or 


^ = Constant 
ds 


and equation (26) can be written 


(30) 


(dt)2 


= 0 


(31) 


which implies that the acceleration tensor 
tions (30) and (22) yield 


p/dy^\ 
dt\dt / 


vanishes along the geodesic. Equa- 


where 

system. 


'dy 


dx 


N+1) 


N+1 _ ( dy 


N+l> 


ds 


O 


ds = xf^ 

ds 


ds = X dt 


(32) 


O 


ds 


(l)| 

O ' MO 


are evaluated at the origin O of the Riemann coordinate 
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reduce to 


The space components of the space -time Christoff el symbols 
the N-space Christoffel symbols, as shown in the remainder of this section. For dt = 0, 


(.%) 


which is the N-space interval. Equation (24) for the space-time geodesic for i 7 ^ N + 1 
becomes, for this case, 


d^y^ . f (X \ dy^ dy^ 

(o- 0 J ds ds 


0 


(34) 


and the equation for the N-space geodesic is 

d^y" , f a dy^ dy^ _ q 

(ds)2 V N 


(35) 


where 


a 

a (p 


N 


is the N-space Christoffel symbol 


and (35) pertain to the same path, an extremal path for 
therefore, 


of the second kind. Equations (34) 
the integral of ds, by equation (33) 


= (3«) 

since dy^/ds is arbitrary at the origin of the Riemann coordinate system, which can be 
anywhere in the Riemann space. 

Raising the first index in equation (28) yields 




’13 


3y3 


+ S' 


1 

k ] 


= 0 


(37) 


and 




i 


17 


Therefore, 




_ 9S^ 
9yi3 


+ S' 




N 


Finally, 


(38) 


(39) 


(40) 



which is equation (27). 


(41) 


The Space -Time Tensorial Derivatives 

A. 

In this section, the covariant and comoving derivatives of a tensor T of any rank 


are expressed in terms of the N-space Christoffel symbols 


a 

jS y 


nate system velocity S^. It is convenient to write the definitions 


and the coordi- 


N 


r . 


,N+1 


1^ 


s 0 


N 



(42) 
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The expression for the covariant derivative of a general tensor is 


8 

®l'--®n|P 8yP 


sr-'®n P/ ®r--®n 


+ . . . + 


m 

q p 


X T 


Sj .-. S ^ 


q '\ rJv • -^m 
Si pj qs2...s^ 


and the comoving derivative is 


q rp 

^n P 


si...s^_iq 


(43) 


ri . . .r 

D rp 1 rn _ rp 


„p 


dt -Si...s^ ‘S1...S IP 


(44) 


By use of equations (36), (38), (39), (41), and (42), equation (43) becomes, for p / N + 1, 


r^ . . .r „ r, . . .r / i i 


yr2 - . .r 


m 


=r--=n|'^ ayP =r--®n \ ' N+iy\^v 


(N+l)r2 ...r / s j r, 

+ T„ ^ 1 - 5„.i “ S + . . . + 1 - 6 


Si...Sn 



N+i; X L" t^N 

=r--=„ =i---=n 


r 
N+1 


){/%} 


m 






- 1-6 


si (s ^ dl is' ■ 'i - is' ' 'i rfsi - • • 

®1 j \®1 Pj ^ ^®2---®n >'®2"-®n ®1 I0N 


( 45 ) 


For p = N + 1, equation (43) becomes 
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The covariant and comoving derivatives of first- and second-rank tensors are suf- 
ficiently important to warrant special expressions. For a contravariant vector 
from equations (45) and (46), 



20 



In particular, 


tN+1 ^ d 

IP ' ayP 


,pN+l 


(49) 


D rpi* _ fTii* i/P — 

-T|pV -T|j,^iV 


= 


|N+1 




9 mr 
— T + 
9t 






dyf- 


'N+1 


y 




l£ 


,N+1 


N 


(50) 


For a covariant vector T„, from equations (45) and (46), 

s ^ 



I T, = T 3 |p vP = V-S = 


X + T 


s|j8 


r/3 


_ 8 _ 

9t 


T - 
s 






- 5 


N+1 1 



T 

y 


^ i t ^ ^ t’lJiN ^ 


(53) 
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r.r„ 

For a contravariant second-rank tensor T from equations (45) and (46), 
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For a covariant second-rank tensor T , from equations (45) and (46), 






N+1' 




5N+1 


J T„ .. -6 


N+1 1 


2 / l ®2 ^ ®2 ^ 






Tc 
N 


( 57 ) 


1 d 


Sj^SglN+l X at Sj^Sg 


(1 - i 


\ "1 j ^ 

L i®iJ 


N 


■ysg 


U _A gN+lNiry^ ^ 


- 6 



nh-1 Ii\- /as^ ^ 


’2 \^/ V 


3^ 




N> 


Sir 


(58) 
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&|0 


r / I* 

For a mixed second-rank tensor Tg [which can be considered to be either Tg or 

T^gV from equations (45) and (46), 



The Rates of Change of the Metric Tensors 


The partial time derivatives of the covariant and contravariant metric tensor com- 
ponents are derived in this section in terms of the coordinate system velocity S^, the ' 


metric tensors, and the N-space Christoffel symbols 



N 
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The covariant derivative of any form of the metric tensor vanishes identically 
(ref. 4, p. 43); hence, by equations (58) and (55), 



Applying equations (12) and (13) yields 



(65a) 


(65b) 


25 



Equations of Motion of a Material Continuum 

In this section, the equations of motion are derived in the Riemann space -time ten- 
sor formalism from fxmdamental principles. This approach, rather than that of merely 
restating the existing equations of motion in the tensor notation, is taken in order to 
demonstrate the versatility and simplicity of the Riemann space-time tensor formulation. 

The equations describing the motions and other properties of a material body 
regarded as continuous are commonly written from either of two viewpoints, the Eulerian 
or the Lagrangian. The Eulerian viewpoint is that of an observer at rest relative to the 
coordinate system; the Lagrar^ian viewpoint is that of an observer at rest relative to the 
material continuum. 

An Eulerian statement of generalized conservation of a general tensor quantity T 
in an inertial rectangular Cartesian system z is 

^(pTV^)+^(pt)=c(pt) (66) 

where p is the mass density and C^pT^ is the per unit volume rate of creation of the 
product of mass M and T, or the source term. Here, "creation" is the increase in a 
quantity from any cause other than transport. The first term on the left side of equa- 
tion (66) is the net per volume flux of MT transported out of a given volume element, 
and the second term on the left side is the observed per volume rate of increase of MT 
within the volume element. 

Equation (66) can be written, by substituting dz^'*’^ = X dt and = X, 

-^(pTV^) = C(pt) (67) 

which generalizes to 

(pfv')|j= C(pf) (68) 

in any allowed coordinate system in Euclidean space -time. The mass continuity equation 
in Cartesian coordinates 

is equation (66) for T = 1 and C(p) = 0; thus. 
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(69) 


(pvO|i = 0 

in Euclidean space-time, by equations (66), (67), and (68). Applying equation (69) to 
equation (68) yields the Lagrangian expression of generalized conservation 

pV't|.=pg(T) = c(pt) (70) 

It follows from equations (A16) and (A17) in appendix A and equation (43) that equation (70) 
holds in any Riemann space-time which is a subspace of Euclidean space-time if all com- 
ponents of and T not contained in the subspace vanish. (All components of 
normal to the subspace vanish, by eq. (A4) in appendix A.) 

Equation (69) also holds in the Riemann subspace if all components of normal 
to the subspace vanish and if 



(71) 


where i is the subspace index, g is the determinant of the covariant spatial metric in 
the subspace, and g is the spatial metric determinant in Euclidean space-time. The 
proof follows from the well-known identity 



and equations (12) for the space -time covariant metric tensor. Let 

2 

® " ■^N+l,N+l ■ ^ 

Then 



■“ , 
1 

^1,N+1 


) 

^1,N+1 



1 

1 

I 

1 

^N,N+1 

+ 

^a/3 ' 

1 

^N,N+1 



^1,N+1 • • • ^N,N+1 

B 


0 ... 0 

1 



(72) 


(73) 
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In the first determinant on the right side, the (N+l)th row is a linear combination of the 
rows; that is, from equations (12), 


gN 

(N+l)th row = — (1st row) + . . . + ^ (Nth row) 

A A 

Thus, the determinant vanishes, and 

A = g (74) 


in any allowed coordinate system in Riemann space-time. Applying equations (74) and 
(72) to equation (69) yields 



= 0 


or, by the assumption of restricted to the subspace. 



(75) 




= 0 


(76) 


which differs from the subspace divergence expression by 



Thus, equation (71) and the assumed restriction of V*" to the subspace are sufficient con- 
ditions for equation (69) to hold in the Riemann subspace. 

Interpretation of the condition expressed by equation (71) is facilitated by letting the 
coordinates y*^, held fixed to define the Riemann subspace, be normal to the subspace 
and, for two y“ coordinates, mutually orthogonal. Then 


S = Sg 


(77) 
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where g is the determinant of the covariant metric tensor in the subspace of the 

coordinates. Substitution of equation (77) into equation (71) yields the condition 
equivalent to equation (71) 



If the subspace of is one -dimensional, then y g is the distance per unit increment 
of y^. If the subspace of y^ is two-dimensional, then '/f is the area per unit cpor- 
dinate increment product dy“^ dy for a = 2,3. 

If equations (69) and (70) both hold in a Riemann space-time, then so does equa- 
tion (68), since 

(pV' f)| . = pV' T|i + (pV')j. f = c(pt) (79) 

Therefore, sufficient conditions for equation (79) to hold in a Riemann space-time are 

(1) All components of T and not contained in the Riemann subspace vanish 

(2) = 0, or equivalently, 

dy^^ ' 'fs ay^ ^ ' 

For some applications of fluid mechanics, sources and sinks for mass are employed 
as a mathematical artifice (see, for example, ref. 5, chapter 5), which allows C(p) to 
be nonzero. No sources or sinks for mass are employed herein, and C(p) vanishes 
everywhere. 

Generalized conservation of mass and momentum is expressible in a single tensor 
equation by setting T = in equation (68): 

(pV^ V^)|. = c(pV^) =f^ (80) 


where, from equation (70), the force per unit volume is 

c(pV“) =f“ 

/ (81) 
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Thus, is a space-time tensor. The Lagrangian form of equation (80) is 


^(pV^) +pV^‘ V|- =f3 (82) 

Either equation (80) or (82) expresses generalized conservation of momentum for 
j 7 ^ N + 1, and mass conservation for j = N + 1. 

The force per unit volume, given by the covariant form of equation (80), 
is separated into the work -producing and velocity -turning components by 



by use of the Eulerian mass conservation equation (69), where is the "space-time 

vorticity" given by 


^ij 



- V 


1 /sV. 



(84) 


Since = 0 (by eqs. (A19) and (A21) in appendix A), an alternate expression for 
is 


Cl). . 


1 ] 






(85) 
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The invariant is, from equation (83), 




( 86 ) 


by equation (70) and the skew-symmetry of cu^j. 

An alternative to the Eulerian momentum equation (80) in terms of kinetic energy 
and vorticity is obtained by partially expanding equation (83) for i = of; 


f = p + 2pV'^ + 2pV^'^^ cu „ 1 

a ^ ay^\2 y ay q',N+1 


= p + 2pV'^ w + 2pX 

9y r 7 ^ ^y® y 


(87) 


where, from equations (18), 


^^N+1 _ „ a Irr 0/8 

3y“ 

Thus, rearranging equation (87) gives 


9V 


pX 


a 


9y 


N+1 


= p — = - P — (I V.\ - 2pV^ co^ + 

9t « ay« \2 7 


9y 


u.s 






(88) 


The Eulerian rate of change of the vorticity in equation (88) is obtained by 

la la 

applying the operators ^ and - to — and — and summing: 

2 9yP 2 9y*^ at at 


1 a _ 1 a ^ _a 


2 9 y /8 \ at / 2 9 vQ! V 9t 


9y" 




a (M _ a fa 1 


^9y 


/3\P 


9y 


a\p 


9y 


3 ‘^ar! * 


9y 


y‘>' oo 
a \ ^Y) 


(89) 
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Thus, if the spatial components of the vorticity field vanish and the acceleration is 
derivable from a potential function, then the spatial components of the vorticity vanish 
for all time. 

The Lagrangian rate of change of vorticity is 


where 





o>, 


ijS \ « j 


_d 

dt 




1 a 

2 ayiS 




+ 










ay 


8V'>' 

9y^ 


+ O) 




8V^ 

9y“ 


by equation (89) and the skew-symmetry of The second bracketed term on the 

right vanishes, by the identity 


Thus, 


and 




= 0 






r 


1 8 
2 ay/3 






(90) 


(91) 


(92) 
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by the skew -symmetry of and the vanishing of by equation (27). The 

source term for the vorticity-mass product 

evaluated by equation (92), can now be substituted into the generalized conservation equa- 
tions (68) and (70), with T = 


CONCLUSIONS 

A space -time tensor formulation for nonrelativistic continuum mechanics has been 
presented which assures that any relationship of quantities expressed as a space-time 
tensor equation is equivalently expressed in all general curvilinear, moving, deforming 
coordinate systems. For example, D'Alambert's "force” is inherently accounted for in 
the tensor expression and transformation of acceleration and force. 

Thus, the space-time tensor formulation applies to situations in which boundaries 
move and deform, and a coordinate surface is placed on the boundary in order to simplify 
the boundary-condition treatment. 

The space -time tensor formulation also applies to situations in which one coordinate 
can be defined as a function of any continuum scalar (e.g., pressure, temperature, or 
density). During numerical integration of the equations ol motion, the continuum scalar 
can then be represented with enhanced resolution by coordinate surfaces which move and 
deform as the scalar varies. Also, the vanishing of the scalar gradient along these coor- 
dinate surfaces can be a significant simplification. 

The generality and compactness of the space-time tensor formulation lend it utility 
as a vehicle for formal manipulations as well. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 8, 1976 


33 



APPENDIX A 


RELATIONSHIP OF CHRISTO FFEL SYMBOLS 
IN RIEMANN SPACE-TIME TO THOSE 
IN EUCLIDEAN SPACE -TIME 

The Riemann space-time considered herein is the locus of constant y®, where 
5 = 3 or 5 = 2,3; and where the y coordinate system is any allowed system over 
Euclidean four -dimensional space-time. Quantities and indices defined with respect to 
the Riemann subspace of Euclidean space-time are denoted by a single bar; for example; 


N = Number of spatial dimensions in the subspace 


O', /3, etc. = 1,. . .,N 



X dt = dy 


N+1 



For arbitrary coordinate displacements in the subspace, the relationships 

Aa,N^-l dy“ dy“ dy« 

%+l,N+l dy^ dy^ 

imply, with equations (12), that 

_ “S 

^5/3 “ ■^5/3 ” ^a/S ” ^5/3 
-^5, N+1 " -^54 > 

%+l,N+l ^ ^44 


(Al) 
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Therefore, from equations (12), 


~ X ^a(p 


O0 _ A_ — = 

^ ■ ^a,N+l 


— ff-— 

X ^OL(^ 


1 

X 




_ — 

a(f) 


or. 



(A2) 


where follows from equation (10), and where indices with a double bar are of 

the y coordinates held fixed to define the subspace. 

From equations (12) and (Al), 


“^44 "" ■^N+1,N+1 ~ (l) ^a/3 ^ 1 


g--s“s^+ 1 


X 


Therefore, 


2g5js5 S^ + gs= S“S^=0 
or, by equation (A2), 

g5jS“S^=0 (A3) 

The left-hand side of equation (A3) is positive -definite; hence 

S“ = 0 (A4) 

Equation (20) expresses the transformation of the Christoff el symbols of the second 
kind between coordinate systems over the same space-time, and so, it is sufficient to con- 
sider the simplifying case in which the y“ coordinates are normal to the spatial subspace 
coordinates and, if N = 1, the y^ coordinate is normal to the y^ coordinate. 

Then, 
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dy3 Ag. dy« 6^ = 0 


dy2 Agi dy“ 6^ = 0 


dy2 Agi dy3 6^3 = 0 


(n = 1 or 2] 

(n = i) 
(n= 1) 


(A5) 


J 


follow from the condition of orthogonality between the general contravariant vectors 
and Q^; that is, E^Q^ = 0. Therefore, 




-^23 “ ® 


(n = 1 or 2 ) 
(n = i) 




(A6) 


and, by equations (12), (A4), and (A6), 


A- - i q:= 
■^Q!4 \ ^Q!0 


i g=-rS^ + ~g=~S^ 
X ^oi(p X ^acj> 


i A — 
A 


= 0 


(A7) 


In general. 




Thus, from equations (A6) and (A7), 


A - A^ - 6^ - A-f A^^ + A== A^^ + A = A^^ 
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or 


si = A|= A^^ 

a PQ! 


(A8) 


If N = 2, then /3 = 5 = 3, and equation (A8) yields Agg ^ 0 and 


a/3^ = o 

A^^ = 0 


(A9) 


If N = 1, then / 3 , a = 2 or 3, and equations (A8) and (A6) yield 
A?^ = 0 




A^^ = 0 


(AlO) 


Thus equations (AlO) hold for N = 1 and 2. 

.a/3 

relationship 


The spatial metric g ^ is uniquely determined by the according to the 


» ^ay y 


(All) 


From equations (12), (Al), and (A6), 


rr-Oi^ ^ r.a/3 - , „ Q! /3 

g S g^jy + g ^ay 


= 8“^* 85.= 6? 

ay y 


= g"^ g-- 
& ^ay 


Therefore, 


-aB aB 
g = g ^ 


(A12) 
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From eqxiations (10), (13), and (A12), 

A®>N+1 = A^* > 

^N+1,N+1 ^ ^44 


(A13) 


The Christoffel symbol of the second kind in the y coordinate system over 
Euclidean space-time (which also holds for Riemann space) is given by (eq. (19)) 



the components of which in the subspace of constant y® are 


(A14) 


By equations (AlO), A^^ vanishes for n = ^; therefore, from equations (A13) and (Al), 
the spatial components of the subspace Christoffel symbol are given by 


where, it should be emphasized, the y^ coordinates (which are those held fixed to define 
the Riemann subspace) are normal to the spatial coordinates in the subspace and, for 
N = 1, normal to each other. Similarly, from equations (AlO), (A13), and (Al), 



(A17a) 

(A17b) 


(A17c) 
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(Al7d) 



(Al7e) 


for the coordinates normal to the spatial coordinates in the Riemann subspace and, 

for N = 1, mutually orthogonal. 

All Christoffel symbols of the form ^ or j j in any allowed coor- 

dinate system vanish, as shown in the following equations, where i, j = 1, . . N + 1. 

In a Cartesian coordinate system z over Euclidean space -time, the covariant derivative 

of 


_ as. 


(A18) 


becomes 


as- 

= — ; = 0 

aSJ 


by equations (18) = 0, = -X in any allowed system^ . Therefore, 



and, by equations (A17), 


(A19) 



Thus, 



(A20) 


(A21) 
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EXAMPLE PROBLEM OF IDEAL FLUID WITH TOP FREE SURFACE 
IN WHICH SURFACES OF CONSTANT PRESSURE 
ARE USED AS COORDINATE SURFACES 

Coordinate System 

Although the governing equations of a material continuum usually have the simplest 
form when expressed in Cartesian coordinates, it is often advantageous in solving the 
equations to set coordinate surfaces coincident with boundaries, in order to simplify the 
boundary treatment. If the boundaries are moving and deforming, then a moving, deform- 
ing coordinate surface is desirable. 

On the other hand, it is advantageous to define coordinate surfaces as surfaces of 
constant pressure, density, 'temperature, or any other continuum scalar variable, and to 
define coordinate lines normal to the surfaces. The coordinate surfaces and lines would 
then move and deform as the continuum variable changes (with respect to a rigid coordinate 
system). The vanishing variable gradient components in the family of coordinate surfaces 
and the orthogonality of the other coordinate lines to the surfaces may simplify the set of 
governing equations. 

These points are illustrated in the example problem of an incompressible inviscid 
fluid with a top free surface moving over the Earth, where, as shown in sketch (b), the 


y^, constant 
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surfaces of constant are the surfaces of constant pressure. Implicit in this repre- 
sentation is the assumption that the pressure increase with depth, which requires only 
that the fluid not accelerate downward at a rate equal to or greater than the acceleration 
of gravity. The y^ coordinate lines follow the pressure gradient and thus are normal 
to the constant -pressure surfaces. The land-fluid boundary is not, in general, a surface 
of constant y^. The land-fluid boundary may rise to intersect the fluid free surface, but 
it is assumed that the slope of the land-fluid boundary is less than the slope of the y^ 
coordinate lines, so that all y^ coordinate lines intersect the free surface. The y^ 
coordinate increments are defined as dy^ = -dp. 

Governing Equations 

For an incompressible, inviscid fluid, the covariant force per unit volume is, from 
equations (12) and (81), 


f = / JP n ^ 

\ dzyj ay« 


9p de 

9y“ “ ^ 9y“ 


(Bl) 


^4(y) 


= A A. . = - f^ g = - f ^ 

(y) i4 X r4> X <P 


(B2) 


where p is the pressure, p is the density, and 9 is the gravitational potential. It 
should be noted that 6 is an invariant scalar in equation (Bl); and, when 9 is considered 
in the inertial z coordinate system, it becomes clear that 9 does not contain the cen- 
trifugal "potential" of Earth rotation which is included in the gravity (as opposed to 
gravitational) "potential." (Neither the centrifugal nor the gravity field is a true potential 
field, because neither is a scalar invariant field.) The centripetal and Coriolis accelera- 
tions appear in the acceleration tensor related to the force per unit volume by 


P 



(B3) 


from equations (70) and (81). 

Equations (Bl) and (B2) in contravariant form are, by equations (13) and (81), 


f 

f 


/3^ 

|3^ 


9P 

.. 

9y 

_9p_ 

9y“ 


- P 

+ P 


d9 \ 



X (p 

g“<^ = 0 




(B4) 
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Expanding equations (B4) by use of 


and 

8yi ay2 


13 23 „ 

g = g = 813 = S 23 = “ 


yields the expressions for 


Uyl ay'^ 


\ay^ dy‘^ i 


f^ = -(-^ + p-^)g^^ 
Vay’^ ay>3/ 


= 0 


(B5) 


(B6) 


where equation (B5) follows from the orthogonality of the y^ coordinate lines to the 


dition 


surface. Equation (82) for j = 4 becomes, by use of the incompressibility con- 
Dp _ 


dt 


0 , 



(B7) 


From equation (28), 


S 



(B8) 


Thus, 




= 0 


or, by equation (49), 

=0 

\y 


(B9) 
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From equations (50)^ (B4), and (B8), 


DU^_DV^_ (1 9p 

dt dt “l^Pay« ay“j 


9t QyY 







= v>' + ^ u>' + f u>'u'^' 

9t ayy ayy V’' 


(BIO) 


Boundary Conditions 

The pressure p at the top surface is the atmospheric pressure p^, which is 
assumed to be constant. The top surface is a y^,y^ coordinate surface; hence 

Q <3 

vanishes and U = S at the top surface. At the land -fluid boundary, the fluid velocity 
relative to the boundary is parallel to the boundary, or 


(u‘-4)"i= (u^-Jb)% = 


(Bll) 


i / i ' 

where jl is the inertial velocity of the boundary 

b dt \ 9t , 




and m is any non- 


vanishing covariant vector normal to the boundary surface. If the region of interest is 
not enclosed by physical boundaries (the top surface and the land-fluid boundary), then 
closure is completed by artificial boundaries, where the boundary condition 


- Jg)n^ = V (B12) 

is applied. Equation (Bl2) expresses the component of fluid velocity normal to the bound- 
ary in terms of the specified scalar v. The slope of the artificial boundaries is 
restricted, as is the slope of the land -fluid boundary, to be less than the slope of the 
y^ coordinate lines; thus, all y^ coordinate lines intersect the top free surface. The 
land-fluid boundary and the artificial boundaries are collectively termed the "bottom” 
boundary, and v is zero at the land-fluid boundary. In summary, the bovmdary conditions 
are 
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and P = Pq 

K - 4 )\ = '' 


(at the top surface) 
(at the bottom boundary) 


> 


(B13) 


Initial Conditions 


Equations (B9), (BlO), and (B13) determine, in principle, the progression of the 
fluid state from the initial fluid state. The initial fluid state can be specified by the initial 
pressure contours in a rigid coordinate system x fixed in the rotating Earth and by the 


initial fluid velocity field U?_\ . For this specification of the initial state, knowledge of 




the density field is not required for integrating the equations of motion, as shown in the 


next section. However, if only the top pressure p^ contour and are given, then 

the density field must be known in order to generate the other pressure contours. A pro- 
cedure for initializing the pressure field for a given top surface configuration, velocity 
field and density field p is discussed in appendix C. |^The vertical velocity 

is determined by the other two velocity components by eq. (B9) and the bottom boundary 
condition. j The y coordinate system is initialized concurrently with the pressure field. 
Subsequently, at the bottom boundary, y^ and y^ are dependent only on x^,x^. 


3 -2 

The x coordinate system is a spherical polar system with x vertical, x 
eastward, and x^ southward; thus, if lunar and solar gravitational perturbations (i.e., 
tidal forces) are neglected, the inertial coordinates z^ can have the origin at the center 
of the Earth, and the velocity of the x coordinate system is 


u 

= J(x) = ° 

(B14) 

h 

= a 

(B15) 


the angular velocity of Earth rotation. 

The bottom boundary vector n^ 

1 2 

boimdary as a function of y and y . 


is determined once for all time at the bottom 


The quantities y^, v, and 
bottom boundary as functions of 



y^ and 


are determined as initial conditions at the 
y^, as are the quantities 
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V 

- ^a(x) 

and, from equation (B12), 

«(y) = •'(y) - ^ • I 


(Bl6a) 


(B16b) 


1 2 

The following quantities, in addition to ^(y)’ ^(y)’ 

initial conditions throughout the region of interest as functions of y^. 


are determined as 
y2, and y^; 


88 ^ 88 dx^ 
8y^ dx^ 8y^ 


(if. = 1,2} (Bl7a) 




(Bl7b) 


U 


3 

(y) 


tt3 





(Bl7c) 


which follows from equation (B9) and the identity 

Q Q 

integration is along y3 coordinate lines and and y^ 
values of U^y) and y3, respectively, 




where the 


are the bottom boundary 





U' 




(Bl7d) 


which is derived in appendix D (eq. (D7)), where the integration is along y^ coordinate 
lines, 


• 5 , 

which is derived in appendix D (eq. (D29)), where the integration is along y3 coordinate 
lines, and at the bottom boundary. 
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Time -Dependent Calculations 


The equations for the time -dependent calculations, equations (B18) to (B27), are 
given in the order of their use in a numerical integration loop. From equations (B5) 
and (BlO), 



^ ^ ^ jjy 

0y2 avV 


ayJ 


9y 


ay 


. ^22 _ 12 _ ^ yV as^ 


ay" 


ayJ 


ay'>' ay 


y 






(B18) 


where v”^ = - S^. During the integration time step, all comoving forces per unit mass 

are held fixed in the z coordinate system. From appendix D (eq. (D33)), 


an 


at^ dy^ Hy) " "^^y^ 




(at the bottom boundary) (B19) 


From appendix D (eqs. (D36) and (D37)), 

.3 


au 


(y) 


at 


rr ^3(y) 
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3 1 ^^3 

^ ! X 


CL 


— lu'^n. + U^n. - V 

at \ 1 2 


(y) 


au 


9y 


m 


- J^) 


(B20a) 


at the bottom boundary, where 



(B20b) 


at the bottom boundary. 


From equations (19), (65), (B5), and (D6), it follows that 
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^^33(y)\ 

I ly* 

haAy)\ 

\ ly’t' 


= ^33 




^ + s'#’ = -v'#’ 

ay3 9y0 



f£33 

3y0 






{a,^,y = 1,2) 




(B21) 


O 

Taking the time derivative of the expression for in equation (Bl7c) yields 



(along coordinate lines) {\p = 1,2) (B22) 


where 


- \ j ^33 (^11 ^22 “ ^12 ^ 12 )- 


From equations (D31) and (Dll) in appendix D, 



(B23) 


12 -12 -3 

since at the bottom boundary, y and y are functions of x and x , and x is a 

1 2 

function of x and x . From equation (D7) in appendix D, 


d q3 

it® 


(y 







(along y3 coordinate lines) 


(B24) 


From equation (D41) in appendix D, 

^ sy 

\ /y0 


(B25) 
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From equation (D29) in appendix D, for = 1,2, 

(along coordinate lines) 



(B26) 


From equation (D43) in appendix D, 


Equation (B27) completes the time -dependent equations. As mentioned before, 
knowledge of the density field is not required for integrating the equations of motion, 
since p does not appear in the time -dependent equations. The density can be recovered 
by means of the three -component momentum equation; from equation (BIO), 


d fde\ 
ysy^j 


33 1 at 9yy 


ay-: 
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ay' 
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y 0 


y TT 0 
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(B28) 


where — ^ is expressed in terms of which are updated by equation (B27) by 

ay ayi dy^ 


Be 

3 


ay 


I 


lP33(y) 

/ BO \ B9 Be 

[bx^J ByP By'P ^ 

^ — 


0,1^ = 1,2) (B29) 


This relation follows from the orthogonality of the y3 coordinates to the y^- and 

2 

y2 coordinates and from g>'0 . (The negative radical in eq. (B29) is 

ayr 8y0 \^ax3 j 

discarded as inconsistent with the condition that the pressure increase with depth.) 
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PRESSURE FIELD INITIALIZATION FOR EXAMPLE PROBLEM 

OF APPENDIX B 


A procedure for initializing the pressure field is discussed in this appendix. The 
velocities u|-y density p(x“); and the configuration of the top free surface 

in the x coordinate system are arbitrarily specified. If the magnitude of the pressure 


gradient were known everywhere, then the configuration of all of the pressure 

contours would be determined by the configuration of the top surface ^of pressure p^), 

since -dy 3/\/^ is the displacement from the contour of pressure p to the contour 
of pressure p - dy^. But, from equations (B6), 


®(y) 


-f 


(y) 


GttO 


9y 


ay- 


ay^ 


(Cl) 


where f^ and p are not known a priori. It will be shown (eq. (C8)) that knowledge 
ay-^ 

of at the top free surface and of the density and inertial velocity fields everywhere 

is sufficient to generate the pressure contours. Then an iterative procedure is suggested 
for determining g^^^ at the top free surface. Once the pressure contours are deter- 
mined, the y coordinate system can be established and the velocity field expressed in 
the y system; then the other parameters can be initialized as discussed in appendix B. 

Taking the space -time divergence of f^ in equations (B6) yields 



ap 89 \ „0Q! _ /„ DU“^ _ 9p DU“ , 

"r * j|o'"9y“ * Ut Jja 


Neglecting the self -gravitation of the fluid so that 



which implies 


m 


N 


g<#>“ = 0, yields 
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,'Pl, 9y® Isy'^' \ > 


i9y^ 
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a 


(C2) 


a 


Substituting “hT” 

P 3y0 dy<P 


( dp \ „cpa _ _j£_ljp^ (l>a [d^\ 

"Mia ' ^\a 


(C3) 


Expandii^ the second term on the right side yields 


dt 


a 


Pi^i. = pK n = K" Hi " ii 


= p U“l K - Sj J = P B Ic U|a 


tQ! 


since ^ ^ ^lo- ^[a = ^ equations (15), (27), (28), and (B9). Thus equa- 

tion (C3) becomes 


/ dp\ = .1 JP__9P_g0» +pu7 uT 

* P3y“9y'#> I'" 1“ 


\ay<p 


a 


(C4) 


Expanding the left-hand side of equation (C4) yields 


_ [JPJ\ = - 


9y 




O' 


9 / 9p \ _ ^ / i 
9y^^ \9y V 9y^ v “ 


g 


(pa 


= _9p. / 3 1 (pa 
" 9y3 \J^ 


3 \ ^(poi 

(p a ^ 


(C5) 


9y 


(/>- > 


4 

(p a 


since = -6^ and / ^ = 0. In equation (C4), 
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£2-g0« = ^g3«_g33 


dy<P 


Hence, 


, ® "\ g'#’“ = 1 j£. „33 . . TjOi 

<!> o,J « p,yS8 +f>U|yO|a 


The covariant derivative of the metric tensor vanishes; thus, 


,3a 


= 0 = ^-^^" (/.) 


N 


= .(g33f!^,^33_L3(if) 


3 

a cr 


9y 


+ g 


o;a / 3 


\[s dy^ 

by use o( g31 , ^32 , ^ {j>'’ ^ ^ j= 


a a 


1 d 


Then 


/g 9y^ 


g . Let 


^12 " ^22 ■ ^12 ^12 


(C6) 


and 


/f=Vi 


33 ^12 


„Q!a 


a a 




2 ag 


33 ^ g33 _1 




/iil 


H- g33 A. ly) 


42 3y‘ 


= -g 


33 I 




(\[^) ^ 33 ^(^ 12 ) 


^12 9y^ 
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Thus, equation (C6) can be written 




1 uf 

^ ay 3 a ^2 ay 3 ^ ay3 ^3 |a |a 


(C7) 


where = v/^ll ^22 ~ ^12 ®12 area per unit coordinate increment product 

1 9 ” 3 3 

dy-*- dy . In terras of the displacement dz = dy along the negative pressure 

gradient, equation (C7) becomes 

K\/^) ^ \l^ dp 

S33PU|^U|^ 


(C8) 


where g^g is expressed in the y coordinate system and pU^^. is an invariant 

scEilar in space-time ^equal to pu|j u||^j . The family of surfaces of constant pressure 

can be approximated from a finite -difference form of equation (C8) by numerical integra- 
tion from the top surface to the bottom boundary, given ggg at the top, the density and 
velocity fields, and the configuration of the top surface. An iterative procedure for 
determining ggg at the top surface is described next. 

The configuration of the top free surface, the velocity components uf-\ and 

2 o 

U^-\, and the density p are specified in an Earth fixed coordinate system x (with x*^ 
as the vertical coordinate). The vertical velocity component is determined by 

u'T = 0 and the bottom boundary condition (u*^ - J^)n_ = v. Then, 


cr 


1 2 

(a) The y and y coordinates are defined at the top surface for the initial time 


in X coordinates. 

(b) At the top surface, ^ggg is initially approximated by 


where dz^ is the displacement element in the y^ direction dy^. This initial 

approximation neglects acceleration a|-'j normal to the top surface and follows from 
equations (B6) with f^ = 0. ^^|z) set equal to zero.j 

(c) The y^, y^, and y^ coordinates and the constant pressure surfaces are con- 
structed in X coordinates from a finite -difference form of equation (C8). The density 
and velocity fields are expressed in y coordinates. 
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(d) The sequence of computations in appendix B consisting of the remainder of the 
initial conditions and the time -dependent operations through equation (B22) for — 
are performed. 


(e) is set equal to a^g^ . The acceleration at the top surface in the 

3 

(z) 

3 


-3 3 

z direction a/ =\ is approximated by (from eq. (BIO)) 


‘(z) - \j^33{y) (“St 


3y'" aya I" 


(y) 


The next approximation to ^ggg at the top surface is, from equations (B6) and (BIO) 
and the relation ggg = l/g^^, 



1 



+ Kp 




where is the previous approximation to ^ggg and K is a constant between 

zero and unity selected by trial and error for fast convergence of ^ggg at the top 
surface. 


(f) Steps (c) to (e) are repeated until a chosen convergence criterion for ^ggg at 
the top surface is satisfied, then the. time -dependent operations of appendix B are 
continued. 
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SUPPLEMENTARY DERIVATIONS FOR TIME- DEPENDENT 
CALCULATIONS OF EXAMPLE PROBLEM 
OF APPENDIX B 


Velocity of y Coordinate System 

AH comoving forces per unit mass are held fixed in the z coordinate system 
during the integration time step; that is (from eqs. (Bl) and (B2)), 


D /I f] = D / 1 ^ 0 

dt\P aj dtl P 9yO! ayo-l 


(Dl) 




' DS'^ _ ^ 

which is automatically satisfied when eq. (Dl) holds, because =01 and 


(D2) 


D/ 1 

dtl P 9yO!, 


(D3) 


D / 90 


dtV ayOf 


(D4) 


Expanding equation (D3) for a = 3 yields, by use of dp = -dy^ and the incompressibility 


condition ^ = 0, 
dt 


1 / 3 

P 13 or 


V^+4f^ + S^ ( . " o)l -0 


Play3 


3 

0 3 


or 


s! - / 3 ^ 

ay 3 ^3 afU 


(D5) 


From equations (36) and (19) and the relation g^^ = 0 for (p 3, equation (D5) becomes 
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^ ^ _ 1 _1_ ^^33 ya 
^ ^33 9y^ 


(D6) 


and the top boundary condition yields, by integration of equation (D6) along 

y^ coordinate lines down from the top surface, 


S^(y\y^y^t) - (y^y^P3,t) 

V / \ O y 2 g33 ay<7 


(D7) 


1 2 

The y coordinate system velocity components S and S are considered next. 
Sketch (c) shows the pressure gradient and pressure contour at a point R of fixed coor- 
dinates y^ at times t and t + At. The directions of the axes z'^ and z^ are 
indicated; however, the point R can move relative to the origin of the z system. The 


H F 




point R is fixed spatially in the y coordinate system. The pressure gradient at 
time t is HG. The z coordinate system is inertial Cartesian with z^ antiparallel 
to HG and the z^,z^ plane parallel to the instantaneous plane of rotation of the pres- 
sure gradient at time t. The pressure contours AB and CD are of the same pres- 
sure (since y3 at R is fixed) at times t + At and t, respectively. The pressure 
gradient at time t + At is FE. The orthogonality of the pressure gradient to the pres- 
sure contours implies that CD and HG rotate about an axis pointing into the plane of 
the paper at the same angular rate co. Therefore, at time t. 



(D8) 


Let be the inertial velocity of the y coordinate system defined in the 
z coordinate system by 
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W(S) = “ 


J 


(D9) 


The definition of is extended to any allowed coordinate system y by requiring 

to be a space-time contravariant vector; thus, 




8z” 



In the y coordinate system, becomes, by equations (9) and (10), 





(DIO) 


(Dll) 


and, by the tensor character of W^, 


^(y) ^(y) 9yi ®(y) ay^ 


wf-x = s? = 0 

(y) (y) ay(7 




(D12) 


From equations (Dll) it can be seen that is the inertial velocity + s\ for 

0 . 

From equation (D8), 


^(y) " 


^ r^ j t ^ r^i =0 
di ^ \ 9t J^a 


(D13) 


or, by equations (D9), 
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'^3|l(2) +^l|3(i) =° 

since the z system is Cartesian. Therefore, 


(D14) 




3azl 


(D15) 


-1 -3 

The z-^,z plane is parallel to the plane of rotation of both the pressure gradient and 
the pressure contour, as specified in the discussion of sketch (c); thus 


9 [9z^ 

9Z^ V 


= 0 


.Of 


9 1 


9z3' 



= 0 




(D16) 


and, by equations (D9), 


Wolo/ = N = 0 


3|2(z) 
''^2|3(z) " ® 


(D17) 


J 


Therefore, 


W.I., ^ 0 

azS 8z2 


w I isl ®y^ - 0 




(D18) 


Combining equations (D18) and (D15) yields 
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Equation (D19) is expressed in terms of the contravariant vector 


by 




-.A-nW^ 


(y) 




(/3=1,2) 


Or, applying 






^=0 

BzP 


(j = 3,4) 


yields 


(A3nW|,^A,„wj'3)^^^g^=0 


(i3,0= 1,2) 


From equations (49) and (D12), 


WT =0 


O 1 o 

Equations (12) and the orthogonality of dy to dy and dy imply 


^3or " ^33 


Thus equation (D21) becomes 






= 0 


(;^,^,0= 1,2) 


(y) 


(D20) 


(D21) 


(D22) 


(D23) 
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which are two equations in the unknowns and 

The determinant of the coefficients dy^/dz^ is the Jacobian for the nonsingular trans- 
formation from y^ coordinates to coordinates, and hence is nonvanishing. Thus, 
the trivial solution is the only solution, and 


^{y) ^%(y) = "^33(y) ^](3{y) = 1,2) (D24) 

or 

= '*'13 (l - «f) 

= -8'^'* 833 

Thus, 

w’jg = ggg ^ 1^2) (D25) 

Expanding the left-hand side of equation (D25) and then the right-hand side yields, by 
equations (19), (B5), and (Dll), 


9W^ 

ayS 


+ 




^33 




= + i w" g w /f!^ . !!3A 

9y2 2 y gyS 9y(/) I 



S33 


ay/3 2 " 8 


®83a 

dy^ 




^33 


8W^ 

dy/3 


_ i 0 ;^^ ^33 

2 ® 
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or 


8 y 3 3 d 


ay/3 


By equations (Dll) and (D26), 


(ft;// = 1,2) (D26) 


833^ 

ay3 gy/3 


At the bottom boundary = y^j, = y^ (x^, x^]; therefore 


-S^= 


4^ 


9t /_0 


= 0 


(ftV/ = 1,2) (D27) 


(D28) 


Hence, 




(;// = 1,2) (D29) 


^yb 

(along y3 coordinate lines) 

Rate of Change of Vector Normal to Bottom Boundary 

The rate of change of the covariant vector n^ normal to the bottom boundary, 
which enters into the application of the bottom boundary conditions, is derived in this sec- 
tion. At the bottom boimdary, 


8n 


0(y) 


at 


'.a 


A/n /-\ 

9tl Qy(p I 


= n_/- 


ct(x) 


zOi 


a fdx^ 


9tUy0 


-Of 

X 


dx‘^\dy^j 9y^\9y*^j 


aV ay'*' 


aV 


iay’/^ dy4> ay4 3y0 dy^j 




1^^ 


dy^ 8y0 \dt J_a ay0 /y/3^ 



(D30) 
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The expression 


( ^z!^ I is the rate of change of the coordinates of a point fixed 

at 

in the spatial x system, or the velocity of the x system relative to the y system 
expressed in y coordinates; thus, by equations (Dll), 


\ / 


(D31) 


Similarly, 


' = W/.-.X - J/,-, 


at (x) ' (x) 


(D32) 


and equation (D30) becomes 


9n 


</>(y) 


9t J_a 


X 




r 


= n 




> ay^ ay^ 






ax^ a 1^ 

■ ay'/' ay'#' 


fy) ■ •’£) 


V(y) ay<t> 


K) - €)) 


(D33) 


Rate of Change of Velocity Component at Bottom Boundary 

From equation (B12) or equations (Bl6), the velocity component U^y^ at the bottom 


boundary is 
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"1 + "2 - ^ 

The invariant scalar n^J°^ is, by equations (B14) and (B15), 


n^J =no/-^S^ 


^2(x) 


Thus, 


9t\ cr ) 


= 0 




Applying eqxiation (D35) to equation (D34) and differentiating yields 

V 

’'l) 

^3 


du: 


M 


at 


-<b 






-<p 


The quantities (-^^1 and ' 


9t - 


X 


4> 




in eqiiation (D36) are evaluated by 


au? A au? V au? \ ^ i 

M = X -( y) = X -- ^ 


9t /.(j, 0x4 


L„cr\ 


ayi dx^ 


M(^\ 


9y^ \ 4® 


9t / (h 


(D34) 


(D35) 


(D36) 


or, by equations (D31), (Dll), and (D28), 
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Rate of Charge of Earth -Fixed x Coordinate System 
Velocity J® in y Coordinates 

The inertial Cartesian coordinate system z, shown in sketch (c), has the origin at 
the center of the Earth; therefore, the velocity of the x coordinate system ^ 

fixed point in coordinates does not change with time; or 


dJ?- 


LEl 


at 


= 0 




Hence, the comoving time derivative of j“ at a point fixed in the y coordinate system 
is 


_ ja ^ dz^ _ ja dz*^ _ ,a ^ ^ / 9z^ \ 

V dt //_N ” |i(z) dt ~ dt ^ l<^(z) I at ) 

N ^(zj \ /yi^ 


or, by equations (D9), 


DJ 


a' 


dt /,.> =j“i(z) "1z) 

(zj 


Thus, by equations (Dll), 


PJ ^ . - .T® . . u/7 . - j. ^ 


a 


dt " |i(y) '^(y) = |a(y) '^(y) = dF(y, 


+ 


0 (J 


(D38) 


or 


f “ J 


rOf 




ay*- 


a 


0 a 




(D39) 


dT® ^ 

by equations (Dll). The term — in equation (D38) is the rate of change of j“ at 


dt 


(y) 


a y grid point, or 
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^*(y) V /y0 

Hence, by equations (D38) to (D40), 

\ /y0 9y^ 

Rate of Change of Gradient of Gravitational Potential 

3 

Since x is vertical, the gradient of the gravitational potential is 


de _ dd ^3 
9xi 9x3 i 


Thus, 


_8_ 

at 





. a / d 0 \ _ - a 


9y 


J3 




dff a^' 

ay'^\ai^ 8y‘^ 



(D40) 


(D41) 


(D42) 


a 

99 

(bA ' 

dy^ 

dx^ 

\ /yOj 


or, by equations (D32) and (D42) and the vanishing of J?-\ , 
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Hence, by equations (Dll), 
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